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;_( ' Abstract 

We study the value of the light quark masses combination niu + rud in QCD using both 
Finite Energy Sum Rules and Laplace Sum Rules. We have performed a detailed analysis 
of both the perturbative QCD and the hadronic parametrization inputs needed in these 
Sum Rules. As main result, we obtain m„ (l GeV^) + m^ (l GeV^) = (12 ± 2.5) MeV 
for the running MS masses. 
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1. The values of the hght quark masses in the Standard Model, are still poorly known. 
We are referring to the masses of the u,d, and s quarks which appear in the QCD Lagrangian 
and which explicitly break the chiral-SU(3) flavour symmetry; the so called "current algebra" 
quark masses. These are the parameters which, among other things, generate the masses of the 
observed low-lying SU(3) octet of pseudoscalar particles i.e.; the Goldstone modes associated 
with the spontaneous breaking of the chiral-SU(3) symmetry. The precise relation between 
quark masses and pseudoscalar masses depends however on the magnitude of the various order 
parameters which govern the dynamics of the spontaneous chiral-SU(3) symmetry breaking. 
The order parameters are genuinely of non-perturbative origin, and this is what causes the 
difficulty in extracting the values of the light quark masses. 

The masses of the vr, K, and 77 particles are however smaller than the spontaneous chiral 
symmetry breaking scale A^, an empirical fact which justifies a priori a perturbative treatment 
of the problem in powers of quark masses. The appropriate field theoretical framework to 
carry out this expansion systematically is called Chiral Perturbation Theory (xPT) Q , Q , Q . 
To lowest order in the expansion. 



(1) 
(2) 
(3) 

Ignoring the 0{mT) contributions, as well as electromagnetic corrections, the expansion to 
lowest order leads to the quark mass ratios estimates [Q]: 
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The modification of these results due to the presence of higher order terms as well as to 
electromagnetic effects, has been an active field of research during the last few years. [For 
a recent review on the subject see ref. [^.] The xPT approach does not fix however the 
scale of the light quark masses. For example, the parameter B which appears in eq. (|l]), 
is proportional to the size of the quark- antiquark vacuum condensate in the chiral limit. It 
is one of the fundamental coupling constants of the effective Lagrangian of QCD at low- 
energies; but its value cannot be determined by symmetry arguments alone. The appearance 
of extra symmetries in the effective Lagrangian ^, preventing the determination of the B- 
like parameters, which appear in xP'^-, within X-P^ itself, reflects the same fact. Ideally, 
the numerical simulations of lattice QCD may fix the values of these B-like parameters to 
a sufficiently reliable accuracy. [For a review of progress in this direction see e.g. ref. Q.] 
Here, we want to re-examine critically the possibility of fixing the value of the running masses 
fnu{s) + rudis) at a specific energy scale s, using general analyticity properties of QCD two- 
point functions. We shall then come back to a discussion of several issues which have recently 
been raised in connection with the determination of the light quark masses. 



2. The relevant two-point function for the determination of the combination of quark 
masses m„ + m^ is the one associated with the divergence of the axial current with the 
quantum numbers of the pion field (see refs. pl to ||l2|,) i.e.; 



$5(g2) = i / d^xe'''-''{0 I T{d^'A'^^-'^\x),d''A^^+'^\0)} \ 0), (6) 

where [q = {u,d,s), and Aj are Gell-Mann's flavour-SU(3) matrices,] 

d^^A(^-'^\x) = (m, + m„)g-(x)i75^i^(z(:E) (7) 

In QCD, the function ^5(9^) obeys a twice subtracted dispersion relation 
(Q^ = —q^ > in our metric: H ): 

^,{q')=^rM+<l'^M+Q'l ^^^^-Im^5(i)- (8) 

The subtraction constants are governed by the low-energy behaviour of QCD. In particular 

^5(0) = -{mu + md){uu + dd) 

^ 2flml{l-5^), (9) 

where 5-,^ denotes higher order corrections. To 0{p^) in xPT Q : 

5^ = Amlf^-^{2Ls - H2) + 0(/), (10) 

to 

where fo denotes the value of the pion coupling to vacuum in the chiral limit, in a normal- 
ization where the corresponding physical coupling is /jr = 92.5 MeV. The constants Lg and 
H2 are two specific coupling constants of the low-energy QCD chiral effective Lagrangian of 
0{p^). The other constant ^5(0) in eq.(P), brings in contributions of 0{p^) so far unknown, 
and therefore will not be used in our analysis. 

At large Euclidean values of the momentum Q^ = —q^, the function ^5 can be evaluated 
using perturbative QCD. The calculation of the first non-trivial order (i.e., two-loops,) was 
made in ref. |^], whereupon the first QCD analysis of the running light quark masses were 
made ( 0, [TO], |jll|, and ||l2|.) The function ^5 is now known at the three-loop level 



in perturbative QCD ( ||13[ and [|lj|). Also, the present determinations of the QCD A- 



MS 



parameter (see ref. |I^ for a recent review,) give values significantly larger than those used 



in the determination made in ref. |12]. These are two of the new ingredients which we shall 
incorporate in the present analysis. 

The basic idea to obtain the value of rriu + md is to use the dispersion relation in eq. 
(|8|) to relate the behaviour of the function ^5 obtained from perturbative QCD, including 
some possible improvement from the knowledge of its leading non-perturbative l/Q^-power 
corrections |16| to specific integrals of the spectral function -Im^5(t). The spectral function 



is an inclusive cross-section which, ideally, could be determined in terms of all the hadronic 
states with overall quantum numbers: I^{J^) = l^(O^), which the operator d^Afj, * (x) 
can produce from the vacuum. The lowest hadronic state contributing to -Im\I'5(t) is the 
pion pole. There is then a gap until one reaches the Svr-threshold. The shape of the low- 
energy end of the Svr-contribution can be determined from lowest order xPT [O] . One would 



like however to incorporate the contribution from higher hadronic states as well. There are 
indeed states, like the 7r(1300) [^, with the quantum numbers of the pion which have been 
observed in spectroscopic studies of the strong interactions. These states contribute also to 
the spectral function -Im^5(t). However, as emphasized by the authors of ref. [0, their 
contribution brings in, a priori, other parameters than those measured in purely hadronic 
reactions. This is yet another issue which we shall also discuss in this new analysis. 

3. A convenient way to encode the information provided by the dispersion relation in 
eq.(P) is to use a system of Finite Energy Sum Rules (FESR's). It has been shown that 
when the free hadronic parameters in a given spectral function are constrained by FESR's, 
the corresponding spectral function automatically satisfies the "heat evolution test" |18| of 
QCD-Hadronic duality. There are three FESR's which are relevant for our purposes: 



f^ at 1 N 

/ --Im*5(i) = -^K(s) + mrf(s)]2s{l + i?o(s)} + ^5(0); (11) 

J^dthm^,{t) = ^[mu{s) + md{s)?ji^l + Ri{s) + 2^^^^y, (12) 

J^^dtthm^i>,{t) = ^[mu{s) + md{s)]^^l^l + R2{s)-^^^^Y (13) 

There is yet another sum rule with two inverse powers of t in the l.h.s. integrand which we 
could consider, but it brings in the unknown constant \I'5(0), and therefore it is not very useful. 
Other sum rules, with quadratic or higher powers of t, become more and more sensitive to 
the behaviour of the hadronic spectral function at high-f values, where there is little reliable 
experimental information at present. 

The right hand sides of the three FESR's above are proportional to Nc, the number 
of QCD-colours. The functions 1 + Rn with n = 0,1,2 denote the corresponding QCD 
perturbative corrections, 

' ""'l-^)D,(t) (14) 
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with D^it) = t^ {^5{t)/t) forn = and 




(15) 
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with \I'5 (t) = -Tpi'^^{t) for n = 1,2, once the sum of the running quark masses: [m^(s) + 
fnd{s)]^, evaluated at the two-loop level and Nc/8tt'^ have been factored out. This D5{t) 
removes the unwanted subtraction constant ^5(0) and ^5 removes both ^5(0) and ^'^{^)- 
Due to the truncated expressions we have for the series of the QCD perturbative corrections 
the explicit form of the i?„-functions depends on the specific choice of the renormalization 
scale |[l^; as well as on whether or not the QCD-running of the two-point ^5-function over 
the circle of radius s in the complex t-plane is taken into account]^. To check the numerical 

*See ref. pc| for a discussion of the relevance of this dependence in the case of the hadronic width of the r. 



influence of the truncation in the QCD series for ^5(t) we shall consider two possible choices of 
the -Rn-functions in the present analysis. One where we resum the leading and next-to-leading 
logs after the integral in eqs. (|l^) and (|l5|) over the circle of radius s is done, by choosing 
the renormalization scale at v'^ = s, and another one where we flrst resum the leading and 
next-to-leading logs by choosing the renormalization scale at z^^ = —t, which we shall call 
"improved," since it takes into account the possible large QCD-running of the two-point \I'5- 
function over the circle of radius s. Other choices like |i^^/s| 7^ 1 in the first case or \iy'^/t\ 7^ 1 
in the second choice do not improve the convergence of the perturbative QCD series for the 
^5-function. 

The quantities (74(04) and Cq{Oq) in the FESR's above, are a short-hand notation for 
the flrst non-perturbative power corrections of dimensions four and six. They are dominated 
respectively by the gluon condensate: 04(04) ~ -^{asC^), and the four-quark vacuum con- 
densate which, in the vacuum saturation approximation ||l^, gives the estimate: Cq{Oq) ~ 
H^ 7r^Qg(gg)^, a result rather similar to the one obtained from the leading behaviour in the 
1/A^c-expansion. In our numerical analysis we have allowed these non-perturbative parameters 
to have values within a rather generous range : 



04(04) = (0.08 ± 0.04)GeV^ (16) 

06(06) = (0.04 ± 0.03)GeV^ (17) 



We use three active quark flavours in the QCD formulae and A-i-/- = 300 ± 150 MeV [15|. 



The procedure we use to extract m„ + rud from the FESR's above, is the same as the one 



discussed in ref. |12]. First, one fixes the choice of the upper end value of s in the hadronic 



integrals, by demanding a good duality between the hadronic ratio of sum rules: 



3 fn dttilm^sft) , , 

^^ad(s)^-^^^^-f— -^, (18) 



and its QCD counterpart: 



l + Rsis) i^^ 


l + R^{s) + 2^^' 
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With s fixed within the duality region, one then solves for rriu + ma using the second sum 



rule in (12); and for (5,^, defined in eq.(lO), from the flrst sum rule in (|r 



4. Another technique to exploit the information encoded in the dispersion relation in 
eq.(^), consists in using sum rules which relate moments of the Laplace transform of the 



hadronic spectral function to their QCD-counterparts. Here, the master equation is |10| 
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where A(M^) denotes the appropriate QCD perturbative corrections, once the overah running 
mass dependence at the M^-scale has been factored out; and the leading and next to leading 
non-perturbative 1/M^-power corrections are explicitly shown. Here again, one has to find 
first a window of duality in the M^-variable. Usually, this is fixed from the comparison of the 
M^-dependence of the ratio of hadronic moments: 



C{M\so) 
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to the M^-dependence of its QCD-counterpart. Then, one solves for the quark masses using 
eg. (pOD , with M"^ fixed in the duality region. Because of the exponential factor, the integral 
in the l.h.s. of eq.(po|) weighs significantly the low energy region of the hadronic spectral 
function. On the other hand, the Laplace transform itself extends from to oo. In practice, 
this means that from a certain SQ-threshold onwards, the hadronic spectral function will be 
identified with its QCD perturbative counterpart. This threshold choice brings in a new 
parameter to be fixed. As shown in ref. |15], the appropriate choice of the onset of the 
perturbative continuum, corresponds in fact to the upper end of the duality region which one 
gets using the FESR's; i.e., sq ~ s. In the present analysis, we shall use the Laplace transform 
technique only as an overall check of the determination of the value of niu + rud obtained from 
the FESR's. Here we have only used the scaling at u'^ = s. The renormalization scale choice 
dependence due to the truncation of the QCD perturbative series is similar to the one found 
in the FESR's analysis. 

5. We shall now discuss the hadronic phenomenological input for the spectral function 
-Im$5(t) which appears in the integrand in the l.h.s. of the FESR's in eqs.([TTD, (|12|), and 

As already said before, the lowest hadronic state contributing to -Im^5(t) is the pion 
pole. There is then a gap until one reaches the Svr-threshold. The shape of the low-energy 
end of the Svr-threshold can be calculated using lowest order xPT, with the result, [A(a, b, c) = 
a'^ + b'^ + c^ - 2ab - 2bc - 2ca]: 
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Notice that the chiral power counting of the pion pole term contribution to -Im^5(t) is 
OljP'); while the contribution from the xPT Svr-continuum is 0{jp^). From the point of view 
of the large-A'c counting rules, the pion pole contribution is leading; i.e., 0{Nc), while the 

^The corresponding expression in eq.(4.14) of ref. |1^ has several misprints which we have corrected here. 



xPT 37r-continuum is suppressed down to 0{1/Nc). Because of the existence of resonance 
states, we expect, however, the three pion continuum to be enhanced, and contribute to 
0{Nc) as weh. The question is how to fold the observed vr' resonance states to the spectral 



function above. The authors of ref. |17] have shown how, ideally, one could measure directly 
the spectral function we are interested in from very high statistics analyses of r — > u^ + Stt 
decays. In the absence of this direct experimental information, it is important to discuss how 
these 7r'(0 ^)-states can contribute, and how to propose an hadronic ansatz for -lui^^it) as 
model independent as possible. 

The situation is in fact rather similar to the one already encountered in the more familiar 
case of the vector-vector correlation function with the quantum numbers of the p. In this 
case, the relevant two-point function obeys a once subtracted dispersion relation: 

U{q^)=U{0) + q' ^- ^-lmU{t), (24) 

Jo t t — q^ -K 

and to lowest order in xPT {0{p'^) in this case,) the spectral function obtained from the 
27r-cut discontinuity is: 

1 1 1 dm^ ■! 

-lmn(t) = —^-{l -—-)-. e{t - Ami). (25) 

Again, this contribution is non-leading from the point of view of the large- A'c counting rules, 
but we know that this is not the full story. The p -resonance also contributes to -lmn(t), 
and its contribution is indeed leading [O(A^c)] in the large-A'^c limit. Knowing the mass My 
and the width Ty (i.e., the /? — >■ 27r coupling gy) from strong interaction experiments, is 
not enough, a priori, to fix the contribution of the same /o-resonance to the spectral function 
ilmn(i). Indeed, at the narrow width approximation, the p-contribution to the spectral 
function: 

ilmn(t) = f^M^7r6{t - M^), (26) 

vr 

brings in the coupling fy of the p to the external vector-isovector current, which is different 

to the coupling gy. These two coupling constants can however be related under the rather 

reasonable assumption of an unsubtracted dispersion relation for the electromagnetic form 

factor of the pion, and p-dominance of the form factor discontinuity. Then, the evaluation 

of the real part of the form factor at q'^ = 0, constrains the two couplings to satisfy the 

relation |21|: 



l = fv9v^- (27) 

The same relation follows, if we parametrize the spectral function -lmn(t) with a Breit- 
Wigner shape function, and we fix the overall normalization to coincide, at the 27r-threshold, 
with the one provided by the lowest order x-fT-prediction: 



-lmn(t) - e{t - 4mJ^^-(l - ^) ^^, _ ^^2 ^ ^2^2^ 



Again, at the narrow width approximation 



-7r5{t-M^), 



and assuming that in the hmit of niassless pions Ty is dominated by the 27r-mode, results in 
the expression 



1 f^ 

-imn(t) = Vttt- 



(30) 



The identification of the two narrow width expressions (|2q ) and ( pOI) leads then to the same 
relation as in eq.(p 



It is well known that the hadronic parametrization in eq.(28) is in fact in very close 
agreement with the low-energy phenomenological spectral function extracted from the e~^e~ 
annihilation cross-section data. Following this analogy, we propose to parametrize the Stt- 
continuum contribution to ^Im^5(t), modifying eq.p^) as follows: 
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The function F which appears as an overall factor in the r.h.s. encodes the modulation due 
to the presence of two vr'-resonancesQ with possible interference (^ is a complex parameter,) 
and with masses Mi 2 and widths ri.2: 
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(33) 



In the hadronic parametrization above we have also allowed for a possible modulation of 
the (1~) 27r-subchannel in the Svr-continuum to take into account the effect of the p(770)- 
resonance. The factor F has been normalized to 1 at the Svr-threshold, but we shall also 
discuss the effect of possibly larger normalization values as suggested in ref. [p^. 



6. The shape of the hadronic continuum contribution to the spectral function in (22) i.e., 
the function 



■''A possible, not yet confirmed, new yr'-resonance has been reported in 
present analysis does not affect any of the results. 



However its inclusion in the 



llm^,{t) ^ e{t - 9ml)-^£^^pf,,,Xty, (34) 

vr HC (16vr^/^)^ 18 

has been plotted in Fig.|l| for s- values between the Svr-threshold and 4 GeV^. The masses and 
widths of the vr'-resonances have been fixed to the estimated values in ref. ||l^ (i.e., in MeV 
units:) 

Ml = 1300 ± 100, Ti = 400 ± 200; Ms = 1770 ± 30, Ta = 310 ± 50. (35) 

The three curves in Fig.|| correspond to the ^-parameter values: ^ = 0.234 + i0.1 (the dotted 
curve); ^ = —0.23 + iO.65 (the continuous curve) and ^ = 0.4 + iO.4 (the dashed curve) 



where the normalization at the Svr-threshold is twice as large as the one in eq. (33). The first 
^-value (dotted curve) is the one which reproduces the best fit to the experimental curves 
which observe the 7r'(1770) in hadronic interactions ||2^. The second ^-value (continuous 
curve) is the one which gives the best duality with QCD. The third ^-value is one which 
gives a good duality with QCD with a normalization at the 37r-threshold equal to 2. Other 
normalizations lead to similar results for the quark masses after imposing the QCD-Hadronic 
duality constraint. 

The duality test is shown in Fig.g. Here, the short dashed curve is the one corresponding 
to the QCD ratio TZqcd{s) in eq.(p9|) and the dash-dotted one is the same ratio evaluated 
with the "improved" set of -Ri_2-functions discussed above. Most of the difference between 
the perturbative and the "improved" version is in R2 and very little in R^. This explains 
why the quark masses obtained show less variation than the duality ratio. The other curves 
are the hadronic ratios, with the same input as in Fig. ^. As can be seen curve 1 has no 
good duality but 2 and 3 are acceptable. We show it further as an illustration of the type of 
variation obtained with the hadronic input. 

Figure shows the results we get for the quark masses. Here we show the value we get for 
the sum of running masses 7n„(lGeV ) -|-mrf(lGeV ) in MeV in the MS'-scheme. The curves 



correspond to the result we obtain solving for rriu + md in the FESR in eq. (p^) , using the "im- 
proved" expression for the QCD-function -Ri(s), and using the three hadronic parametriza- 
tions shown in Fig.|^. The continuous curve corresponds to the hadronic parametrization 
which passed best the QCD-Hadronic duality test. For this case we also show the result with 
the perturbative expression for Ri. The stability of the result is rather remarkable. If the 
same good QCD-Hadronic duality is wanted, then a modification of the hadronic parametriza- 
tion of the Svr-continuum in (^2|), including, for instance, a global normalization factor at the 
37r-threshold to the two vr'-resonances overall function F in (p3) (as suggested in ref. [|l^]), or 
varying the complex parameter ^, leads to a spectral function very similar to the continuous 
curve in Fig. ffl and therefore to the same results from the FESR's. As an example we have 
shown the effect of a normalization factor of 2 and a variation of ^. Thus, we quote the results 
for the best duality hadronic parametrization: 

m„(lGeV2) + md(lGeV2) = (12 ± 2.5) MeV, (36) 

where the error reflects changes in the input parameters. 

7. It is well known that the scalar channels can be affected by instanton effects, however, 
there is no clear computational scheme to get a reliable estimate of these effects (for an 



attempt see [^). Nevertheless, it is also obvious from the results in this reference that if the 
duality region is for energies around s ~ (2.5 ~ 3.5) GeV^ as we have, these effects can be 
absorbed in the quoted error bars for the masses above from the FESR analysis. 

Figure |^ shows the result we get for the parameter (5,^ defined in eq. (|9|) . This is the result of 
solving for (5,^ in the first FESR in ([TTI). Again the continuous curve is the one corresponding 
to the hadronic parametrization which passed best the QCD-Hadronic duality test i.e., the 
one which leads to the continuous lines also in Figs.|^ and y. For the QCD function Rq we 
use the perturbative expression. The improved one is within the errors. The stability of the 
result is less impressive than for the quark masses in Fig. y, but good enough (notice the 
vertical scale,) to obtain an estimate of this parameter. 

6^ = (3.5 ± 1)%. (37) 

From this result and using ( [l0| ) we can obtain for the scale independent quantity [^ 2L8 — H2 

2L8 - H2 = (2.9 ± 1.0) X 10-3 , (38) 

and from (|9|) 

{uu + Jd)(lGeV^) ^ _^^^^3 ^ o.003)GeV3. (39) 



This result is for the non-normal ordered quark condensate (see refs. [25] and p 



We can also combine our result for the light quark masses in ( |36|) with the recent results 



for the strange quark mass in refs. |25, 26| and |M| to get the ratio of quark masses (we use 
m,(l GeV^) = (175 ± 25) MeV): 

TTI 

r = ^ = 29±7. (40) 

m 

Although the error bars are large, this result is free of the uncertainties noted in Q since it 
is obtained from QCD. 

We can also combine our result with the estimation of the next-to-leading corrections to 
Dashen's theorem in refs. [E^ and VOl 



{Ml, - Mio)^^ = (1.9 ± 0.4) (rnl, - ml^)^^ . (41) 

This result translates into, using the ratio from [^, 

m^-m^ ^ rn| (Mio-iV4)Q,p^2 ^2 ^ ^ ^ ^^_^ _ 

nid + ruu Ml Ml - ml Am^ V / ^ ^ 

and with the values for the quark mass ratio r above we get 

TTI 

—^ = 0.44 ± 0.22. (43) 

rrid 

We shall next discuss the consistency of our results for the quark masses, using FESR's, 
with the results from the Laplace transform technique. With the replacement in eq. ( pll) 
incorporated in the spectral function in (|2^), the master equation in (pO|) can be written as 
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(44) 



with P{M'^, s) the corresponding perturbative QCD expression given in the appendix. Solving 
for the quark masses in this expression results in the plots shown in Fig. |^. As in the previous 
figures we show the three hadronic parametrizations of Fig. || for s = 2.5 GeV . The results 
for s = 3.5 GeV are very similar. The agreement between the two determinations is very 
good. We have also checked that the hadronic parametrization gives a reasonable duality for 
the ratio in eq. (|2T|). 



Acknowledgements 

EdeR thanks NORDITA where part of this work was done for hospitality. JP thanks the Leon 
Rosenfeld foundation (K0benhavns Universitet) for support and CICYT (Spain) for partial 
support under Grant No AEN-93-0234. This work was also done under NorFA grant No. 
93.15.078/00 and INTAS 93-283. 

Appendix 



1. The Two-Point Function ^5(9^) in perturbation theory. 



With z^^ the renormalization-scale, and Q^ = —q^ 
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The non-trivial coefficients are 60 = Ij &i = 17/3 and 

'65 2 



10801 _ 39 

144 2 ^^ ^ 



24 



C(3) nj, 



(A.2) 



with Uj the number of active quark flavours. The other coefficients are fixed by renormal- 
ization group properties with the results: c\ = —71/2; c^ = (6i/4)(/3i — 271) — 72/2 and 
t^2 = (— 7i/12)(/9i —271), where /3j and 7^ are the coefficients of the perturbation theory series 
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expansion in powers of agiy^^JT^ of the /?(as)-function associated with the QCD couphng con- 
stant renormahzation: /3i = — ll/2 + nj/3 and (i^ = — 51/4 + 19n//12; and the 7(as)-function 
associate with mass renormahzation: 71 = 2 and 72 = 101/12 — 5nj/18. From eq. ( [A.l| ) we 

can easily obtain the expressions for -D5(t) and ^5 (t) with t = q"^, needed in eqs. ( |T^ ) and 
([T5|) for the evaluation of the FESR's. 



2. The Spectral Function -Im^5(t) in perturbation theory. 



-Im*5(i) 

TT 



QCD 



^K(.2)+m,(.2)ft|l + ^ 



61 + 2ci log -^ 



+ 



«s(^') 



vr 



62 + 2C2 log ^ + ^2(3 log^ -^ - vr^) 



+ o 



(A.3) 



3. The running coupling constant and the running mass. 



With ai = 2/ f — /?! log ^j , the running coupling at the two-loop level, at the scale Q^, is 



a{Q ) = — - — = ail 1-ai— loglogp- 



(A.4) 



At that level, the renormahzation group summed expression for the running mass at the 



Q -scale is: 



m{Q ) = m a{Q 



-h 



1 + f a(Q" 

Pi 






(A.5) 



where m fixes the arbitrary constant of motion of the renormahzation flow differential equa- 
tion, which defines an invariant mass. 

4. The Laplace Transform M^{M'^) in Perturbation Theory. 



The operator L which transforms \&5 (Q^) into the function 



M^{M^ 



{M^Y JO 



00 _ t I 

cite M^ -Im^5(t), 

vr 



IS 



L = lim 



2\N 



N^oo 1 r(A^) (9Q2) 



with M2 = Q'^/N finite. 



(A.6) 

(A.7) 
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As explained in the text, in practice we have identified the hadronic spectral function 
-Im^5(t) in the r.h.s. of eq. ( A. 6 ) with its perturbative QCD expression from some threshold 



sg — s onwards. This allows us to rewrite eq. ( [20| ) as 



(M^)'^ J ^^^ ^ -Im^5(i) = -^[mu{s) + md{s)Y 



with P(M2, s) defined by 

1 f , - * 1. 



/•S J 1 AT f ^ 

/ dte-^ -Im^5(t) = -^[mu{s) + md{s)f\l + P{M\s)\. (A. 

Jo TT QCD StT^ I^ J 



(M2) 
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Figures 



Hadronic Spectral Function 




s [GeV"2] 



Figure 1: Shape of the hadronic continuum contribution to the spectral function in eq. ( P2[) 
in units of 10~^ GeV . The continuous curve is the one which gives the best duaUty with 
QCD. The dotted curve is the one which reproduces the best fit to the experimental curves 
which observe the 7r'(1770) in hadronic interactions. The dashed curve is the one with the 
best duality when the normalization at the Svr-threshold is twice as large as in eq. (p^)- 
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Figure 2: QCD-Hadronic duality test. The short-dashed curve is the one corresponding to 
the QCD ratio TZqcd{s) in eq. (^) with u'^ = s. The dash-dotted curve is the "improved" 
ratio. The other curves are the ones corresponding to the hadronic ratio TZhad{s) in eq. ( |l8|) 
using the same input values as in Fig. |^. 



Light Quarl^ Masses 




s [GeV**2] 



Figure 3: The running quark masses mu{lGeV ) +mii{lGeV ) in the M5-scheme. The three 
curves correspond to the results we obtain solving for rriu + rud in the FESR in eq. (|l2[), 



using the "improved" expression for the QCD-function Ri(s), and using the three hadronic 
parametrizations shown in Fig.|^. The short-dashed curve shows the change if we use for 
the hadronic input the 2nd parametrization and the perturbative expression for the QCD 
counterpart. 
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Figure 4: The parameter 5-,^ defined in eq.(^) which results from solving for 5-,^ in the first 
FESR in eq.(|ll|). The three curves correspond to the hadronic parametrizations shown in 
Figjl| and using the perturbative QCD expression for Rq. 



Mass from Laplace 
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Figure 5: The running quark masses 7n„(lGeV ) + r7T,rf(lGeV ) in the MS'-scheme, obtained 
with the Laplace transform sum rule in eq.(H3). Curves 1,2 and 3 correspond to the hadronic 



parametrizations shown in Fig.|l| and for s = 2.5 GeV . The results for s = 3.5 GeV are very 
similar. 
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